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Abstract Collective-density variables have proved to be a useful tool in the 
prediction and manipulation of how spatial patterns form in the classical 
many-body problem. Previous work has employed properties of collective- 
density variables along with a robust numerical optimization technique to 
find the classical ground states of many-particle systems subject to radial 
pair potentials in one, two and three dimensions. That work led to the iden- 
tification of ordered and disordered classical ground states. In this paper, we 
extend these collective-coordinate studies by investigating the ground states 
of directional pair potentials in two dimensions. Our study focuses on di- 
rectional potentials whose Fourier representations are non-zero on compact 
sets that are symmetric with respect to the origin and zero everywhere else. 
We choose to focus on one representative set which has exotic ground-state 
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properties: two circles whose centers are separated by some fixed distance. 
We obtain ground states for this "two-circle" potential that display large 
void regions in the disordered regime. As more degrees of freedom are con- 
strained the ground states exhibit a collapse of dimensionality characterized 
by the emergence of filamentary structures and linear chains. This collapse 
of dimensionality has not been observed before in related studies. 

Keywords First keyword • Second keyword • More 



1 Introduction 

The Fourier components of a set of position or density variables are commonly 
referred to as collective-density variables and the resultant space is called the 
collective-coordinate space. The sort of analysis involving collective-density 
variables naturally arises in the context of the A'^-body problem, where a 
collective-coordinate transformation will often simplify seemingly intractable 
problems. The reader is referred to Feynman's discussion of the superfluid 
properties of ^He or a series of papers of Percus and Yevick for illustrative 
examples of their use [HOIS]. 

Recently, there has been renewed interest in collective-density variables 
and the associated ground states for a class of soft, long-ranged pair poten- 
tials [4l[5ll6l[7l[8ll9l ll0 l[T2l[TT| . A proof of the existence of periodic ground states 
for this class of soft, long-ranged, oscillatory pair potentials utilized such an 
approach 6 , 7 . Duality relations were derived, connecting these ground states 
to those of particle systems subject to certain short-ranged, soft potentials, 
allowing a rigorous understanding of the ground states of such many-particle 
systems PIT^. Various studies utilized properties of collective-density vari- 
ables to demonstrate the existence of both disordered and crystalline ground- 
state configurations for classical many-particle systems subject to isotropic 
pair potentials in one, two and three dimensions [4l[5ll8l [T0llll) . These interac- 
tions also exhibit unusual excited states; see Ref. [IH] and references therein. 

The present paper extends these last collective-coordinate investigations 
to such classical systems in subject to anisotropic or directional potentials. 
As shall be seen, the introduction of anisotropy results in several striking, 
unintuitive features of the ground state configurations. In the disordered 
regime, ground states display a tunable void structure. Ground states exhibit 
a collapse of dimensionality characterized by the emergence of filamentary 
structures as more collective density variables are minimized. 

It is useful to remark on the possible applications of the results presented 
here. The pair potentials we study are directional in addition to being long- 
ranged, oscillatory, and bounded at the origin, characteristics which are in 
some sense artificial. However, such interactions might present themselves in 
polymer systems J-4,15,16 or in systems of particles adsorbed on a metal- 
lic substrate in which the conduction electrons directionally interact with 
the adsorbed material |17) . Furthermore, a decoration of disordered classi- 
cal ground-state configurations generated from the collective-coordinate ap- 
proach has found real- world applications in the design of disordered dielectric 
metamaterials with large and complete photonic band gaps [18\ . 
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In the present paper we use "directional interaction" to refer to pair 
potentials between point particles that depend on direction relative to a 
fixed spatial coordinate system. Thus, we do not consider cases of anisotropic 
particles with rotational degrees of freedom, such as "patchy colloids" [HI 
[20]. 



2 Definitions and Preliminaries 

We proceed to define the collective-density variables of a particle distribution 
and elucidate some of their important properties. We are given a configura- 
tion of N point particles confined to a volume f2 in d-dimensional Euclidean 
space M'', subject to periodic boundary conditions. The local particle density 
of the system can be expressed as a sum of delta functions: 

N 

p(r^) = ^<5(r-r,). (1) 
The collective-density variables are then conventionally defined as: 

N 

p(k) = ^exp(ik.r,). (2) 

which, to within a constant factor, are the Fourier components of Eq. (1), 
where Vj denotes the j"^ particle position and k is a wave vector defined by 
the parameters of the d-dimensional fundamental cell. For a rectangular cell: 

/ 277711 27rn2 2Trnd\ 
k= — — ,— — — , m, ... e /. (3) 

\ ^1 1^2 / 

In principle, there are an infinite number of p(k) at our disposal. However, 
since the system only possesses dN configurational degrees of freedom, the 
point pattern is completely described by dN collective-density variables (al- 
though it has been shown that as few as half as many may completely con- 
strain a system, at least in one dimension [T). 

There are some properties of the p(k) , stated now, which will be useful 
in the following analysis. For the zero k- vector, the corresponding collective- 
density variable takes on a value independent of particle positions, which is 
simple to compute from the definition above: 

m = N. (4) 

For all other p(k), we have the following properties which follow from the 
relations described above: 



p(-k)=p*(k) 
0< |/5(k)| < TV. 



(5) 
(6) 
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As will be made clear in the following development, it will be useful to define 
the "real collective-density variables" in the following way: 

C« = (7) 

N 

= cos[k.(r,-rO]. (8) 

The real collective-density variables have the following properties: 

C(0) . ^ (9) 
C(-k) = C(k) (10) 

4 < cw < (u) 

We can also relate the real collective-density variables to the structure 
factor S'(k) which describes how a particle configuration scatters incident 
radiation. In terms of the real collective-density variables, 5'(k) is shown to 
be 0: 

We now introduce w(r), a pair potential which is symmetric under inver- 
sion and translationally invariant. In addition, we assume that its Fourier 
transform, y(k) is well-defined, such that: 

V{k) = f dr v{r)e^-'' (13) 

vir) = ^Y.^{k)e-^''''. (14) 

k 

Then, the total interaction energy of the system is given by the sum of the 
interaction potential over all particle pairs, which we can relate to the C(k): 

N 



^>(r~)= ^ «(r, -rO (15) 

^^E^(k)[p(k)p(-k)-iV] (16) 

k 

= ^Y.^ik)C{k). (17) 



k 



We proceed to choose an interaction whose Fourier transform is strictly 
positive on a finite set of wave vectors. A, and zero everywhere else. We 
then observe that by minimizing the C(k) for all k S A, we will obtain an 
absolute minimum for the total energy of the system. As shown earlier, the 
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C(k) are inversion symmetric, so constraining C(k) to some value necessar- 
ily constrains C(— k) to that same value. In constructing our optimization 
procedure, then, it is only necessary to consider half of k-space and to de- 
fine our potential on one half of k-space. We can then relate the number of 
constrained C(k) to the total degrees of freedom of the system, defining a 
dimensionless quantity [5j: 

- ^ 

where M is the total number of independent C(k) constrained, d is the di- 
mensionality of the particle container, and TV is the number of particles. 



3 Computational methods 

In all of the following simulations, we examine particle configurations in 
confined to square fundamental cells, the dimensions of which have been nor- 
malized so that the total particle density is unity. Periodic boundary condi- 
tions are imposed on the system. In order to remain consistent with previous 
work, we try to limit our simulations to iV = 418 particles, when appropri- 
ate. This number was chosen for preliminary studies since 418 particles can 
be arranged into a nearly undeformed triangular lattice when confined to a 
square fundamental cell [5]. 

In most cases, the initial conditions for the simulations are Poisson dis- 
tributed configurations drawn from the same distribution. For these instances, 
the simulation procedure is analogous to an "infinite quench" in which the 
system is brought from an essentially infinite temperature (and a maximally 
random configuration) to zero temperature (and an energy-minimizing con- 
figuration), instantaneously. For x > 0.5, convergence to an absolute mini- 
mum from Poisson initial conditions becomes much more difficult (with only 
about one in twenty runs converging in the best cases). In order to ensure 
convergence to a global minimum, randomly perturbed lattice patterns cho- 
sen by ansatz were used as initial conditions. 

The minimization itself was accomplished with an implementation of the 
MINOP algorithm ^21j . It was shown in previous work that MINOP is advan- 
tageous when compared with other optimization schemes [10]. The conjugate 
gradient method was used in preliminary trials, but was later shown to be 
much less likely to find a global minimum, although the procedure would 
terminate in marginally fewer steps using this method |5]. 

Because of errors resulting from the numerical imprecision of the com- 
puter simulations, ground-state energies were never found to be exactly min- 
imal in practice. The cutoff for determining whether or not a configuration 
was a ground state or not was set to be 10"^*^ natural units from the minimum 
[which can be calculated from Eqs. (11) and (17) and the form of the po- 
tential], although energies were predominantly much closer to the minimum 
than this (on the order of 10~^° natural units). 
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4 Directional Pair Potentials 



Previous work has characterized several results for isotropic potentials in 
two and three dimensions [5,8,10 , which correspond to finite, hyperspherical 
regions in k-space on which ^(k) is a positive constant (specifically set to 
unity in our simulations in order to simplify analysis), i.e.: 

^^^> -jo, K< |k|, ^^^> 

where K is some positive constant. Any ground-state configuration for this 
isotropic pair interaction has been termed "stealthy" [TU] because the asso- 
ciated structure factor S{k) is exactly zero for k < K, i.e., single scattering 
of incident radiation is completely suppressed for k < K. Moreover, by con- 
struction, such stealthy configurations are also "hyperuniform" [22], since 
infinite-wavelength density fluctuations vanish. In what follows, we will refer 
to this region on which ^(k) is non-zero as the "k-space exclusion region," 
since C(k) (or, equivalently, S'(k)) is minimized on these wave vectors. In 
the infinite-volume limit, as k-space becomes dense, the real space potential 
becomes: 

vir) = ^MKr), (19) 

where Ji is the first order Bessel function of the first kind and r denotes the 
radial distance from the particle. This potential is both finite at r = and 
long-ranged and oscillatory and as r — >■ oo: 

Utilizing this class of potentials, Uche et al. 5 categorized the emer- 
gence of three ground-state regimes in two dimensions as a result of in- 
creasing X from zero: disordered, "wavy crystalline," and crystalline (dis- 
ordered and crystalline configurations were discovered in three dimensions, 
but wavy crystalline configurations were not, though they are postulated to 
exist). Disordered ground states were found to minimize the potential en- 
ergy for X < 0-5 and are novel for classical systems (although disordered 
ground states for quantum mechanical systems have been shown to exist 
[33]). Wavy crystalline ground states appear for the range 0.5 < x < 0.7 (al- 
though this was only demonstrated for two-dimensional systems). Perfectly 
ordered crystalline states appear as the k-space exclusion region comes close 
to the boundary of the second Brillouin zone, for x 0-7 in two dimensions. 

It was found that if x is large enough, the sole ground state for this class 
of potentials in two dimensions is the triangular lattice. Initial theoretical 
studies proved this to be the case [SMZIEIEI] • Later numerical studies utiliz- 
ing the collective-coordinate approach confirmed the theory. Namely, it was 
shown that above a certain critical X = Xc, the unique ground state of a 
potential is the triangular lattice [8lll0|. Below this critical value, however, 
the ground states remain degenerate. 
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4.1 Sample directional potentials 

In the present paper, we analyze the ground states for several anisotropic or 
directional potentials with various k-space exclusion regions. For example, 
wc studied the case where the exclusion region is an ellipse. Wc explored 
the parameter space of the model varying both x^ defined in the previous 
section, as well as b/a, the aspect ratio of the ellipse. Various quantities of 
import were calculated, including the static structure factor, S'(k), which, as 
previously stated, describes the intensity of scattered waves, as well as the 
pair correlation function, .92(1"), which describes the conditional probability 
of finding a particle at position r given a particle fixed at the origin. The 
numerical and analytical results for the ellipse potential are presented in the 
Appendix. 

Two other potentials in two dimensions were analyzed, which will be men- 
tioned here, briefly: the k-space square and the k-space plus sign. Disordered 

ground state configurations for these potentials did not possess any notice- 
ably unique properties. Ordered ground state configurations were not found 
for either potcintial. Thcisc; potentials were not pursued any further, because 
their ground states did not display any qualitatively new features. 



4.2 Two-circle directional potential 

We focus our attention on the analysis of the ground states for directional 
potentials whose k-space exclusion region is two circles symmetric about the 
origin. We define the dimensionless parameter, 



IT' 



(21) 



the ratio of the circles" displacement from the origin, k^isp, to their radius, 
K. Once again, the parameter space was explored, varying both x and a. 
The quantities <S'(k) and 52(1") were also calculated. See Figs. 1 and 2 for 
visualizations of this potential. 



4.3 Simulation results 

Several ground state configurations were generated for all potentials, over 
a variety of parameters. As was the case in previous studies, convergence 
for X < 0.5 was trivial for all potentials. For the two-circle potential, on 
the other hand, convergence was obtained for all x up to x ~ 0-74 and 
for the entire tested range of the parameter a. Three difli'erent regimes were 
identified: disordered (with void regions), "filamentary," and "ordered" (with 
a collapse of dimensionality); see Fig. 3 for examples of such configiirations. 
In general, it was observed that convergence was slower for the two-circle 
potential, likely having to do with the "roughness" of the energy landscape, 
as discussed in the next section. Sample ground-state configurations for the 
two-circle potential, ranging from disordered to ordered, are presented below. 
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Fig. 1 Two-circle potential, Fourier (left panel) and real-space (right panel) rep- 
resentations. The dimensionless ratio of displacement distance to circle radius, 
a = kdisp/K, is 1.4548. The fraction of constrained degrees of freedom, X) is set 
to be 0.4461. For these values, the particle soft core has an a;-direction diameter 
of approximately 1/20 the box length and a y-direction diameter of approximately 
1/10 the box length. The energy is scaled such that Vo in Eq. (26) is unity (Color 
figure online) 



8 = 9 = 7i/4 e = 7i/2 




5 10 5 10 5 10 



r r r 

Fig. 2 Real-space two-circle potential Eq. (26), v{r,9) plotted radially (vs. r) at 
fixed angle, 6. The dimensionless ratio of displacement distance to circle radius, 
a = kdisp/K, is 1.4548. The fraction of constrained degrees of freedom, X) is set to 
0.4461. The difference in the x and y-direction soft core diameters can be observed. 
The energy is scaled as in Fig. 1. 



5 Two-circle potential: translation of constrained C(k) 

We go on to derive relations between the two-circle potential and the isotropic 
"stealthy" potential in two-dimensions, which will aid in the analysis of our 
results. In two-dimensional k-space, we define the indicator function for some 
compact set K. centered at the origin to be: 



(22) 
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Fig. 3 (a) Disordered, (b) "filamentary," and (c) ordered configurations for tlie 
two-circle potential (Eq. 26), with fixed kdisp = 4.8710 (in units of inverse length) 
and a = 1.4548 (x = 0.4461), a = 1.3439 (x = 0.5227), and a = 1.1384 (x = 
0.7285), respectively, with A*' — 418 particles. Within each row of the figure, the left 
image is a sample configuration and the right image is the logarithmically plotted 
structure factor ^(k). Black regions indicate that 5(k) has been constrained to a 
minimum. Lighter regions indicate that S(k) > 1. (Color figure online) 



We assume that this indicator function has a well-defined Fourier transform 
/x(i")- We proceed to translate the indicator function by some constant vector 
kdisp and its additive inverse — k^isp, giving us 0k{^ — '^disp) and 0kO^ + 
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kdzsp), respectively. We can signify the union of the translated sets by the 
sum of the indicator functions. By the linearity of the Fourier transform 
operation, the Fourier transform of the union is then the sum of the Fourier 
transforms of the two translated indicator functions. The shift property of 
the Fourier transform [2 5) gives us: 

J-r[/(k - krf,,p)] = J-r[/(k)]e-2-'^-- "- (23) 
Putting this together, we see: 

= 2cos(27rkd,,p •r)/^(r). (24) 

Therefore, for the two-circle k-space representation, the real-space po- 
tential will be that of the isotropic pair potential, modulated by a factor of 
2cos(27rkdisp • r): 

v{r) = -2 — Ji{Kr)cos{2TT'kdtsp ■ r) (25) 
zvrr 

= Ji{Kr)cosi27raKkd,,p-r), (26) 

where we have switched from radial to vectorial representation and where 
kdisp is the unit vector in the direction of kdisp- The parameter K — a/x/tt 
is the radius of the transformed circle. The parameter a is as defined in Eq. 
(21) and the other parameters are as defined in Eq. (19). 

We observe that the pair correlation function is equivalent to that of the 
ground states of the isotropic stealth potential for fixed angle 6 = tt/2. As 
the angle goes to 6* = 0, we see the development of a large peak near r = 
and the emergence of long-ranged oscillations; see Fig. 4 for an illustration of 
this behavior. It was mentioned previously that the optimization procedure 
took significantly more time to find a ground state for all x for the two-circle 
potential. It can be hypothesized that this is the case because the added 
modulation makes the pair potential as well as the total energy landscape 
more "rough" (i.e, many local minima), thus, making the search for the 
absolute minimum more difficult. 



6 Discussion 



6.1 Disordered ground states 

Ground-state configurations were generated at fixed x values (x < 0.5) while 
varying the parameter a for the two-circle potential. For a > 1, it is apparent 
that the pair correlation function is modulated by a factor of cosine in the x- 
direction (see Figs. 6 and 7 for pair correlation and structure factor plots). In 
addition, if x is chosen from the approximate range 0.3 < x < 0.5, noticeable 
void regions begin to emerge (calling to mind, though qualitatively, the voids 
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e = e = 7t/4 e = 7t/2 




r r r 



Fig. 4 Plot of the pair correlation function 32 {r, 9) versus r for selected angles 
{6 = 0, 7r/4, 7r/2), x = 0.4461, a = 1.4548, which has been averaged over 1000 con- 
figurations. For 6^ = 0, the pair correlation function displays long-ranged, oscillatory 
behavior which disappears as 9 — > ■k/2. 



present in the galaxy distribution |26jV It was found that the shape and size 
of these void regions could be tuned by varying both parameters of the model. 
For a ^ 1, the void regions display defined, polymer-like boundaries, which 
seem to be directed along the a:-axis. As a is increased, the void size stays 
approximately the same, while the distribution of particles around the voids 
seems to become more uniform (see Fig. 5). This seems to be in accordance 
with the nature of the potential. For smaller a, the cosine modulation creates 
deep potential wells along the cc-direction (refer to Fig. 2 for a visualization 
of this). As a is increased, the oscillations of the cosine modulation become 
much more rapid along the axis, decreasing the width of the deepest wells, 
while introducing smaller wells with greater frequency in the x-direction. It 
is useful to note that for a > 1, the two-circle potential is marginally stable 
(its volume integral or, equivalently. V{k = 0), is zero [27]), which may help 
to explain the existence of void regions. However, this marginal stability 
property is a necessary but not sufficient condition for the appearance of 
voids, since they only arise for sufhciently large values of x when a > 1. 



6.2 Ordered ground states 

The ordered ground states of the two-circle potential system display surpris- 
ing characteristics. For 0.5 < x < 0.7, ground-state configurations exhibit 
an increasing collapse of dimensionality with increasing x- For x near 0.5, 
the particles fall into polymeric chains which are oriented in the direction 
of displacement of the circles in the k-space representation of the potential. 
When X is increased to its maximal value (x ~ 0.7), the particles collapse to 
several nearly perfectly linear chains, once again oriented in the direction of 
displacement. The chains are randomly distributed along the y-axis and the 
number of linear chains increases as the parameter a is increased. In some 
cases, the chains take on qualitatively sinusoidal shapes which are ordered 
along the x-coordinate. 

We now clarify exactly what we mean by "collapse of dimensionality," re- 
ferring to Fig. 3c. Each one-dimensional linear chain (along the x-direction) 
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a = 1 .9397 



a = 2.9095 



a = 3.8794 



Fig. 5 Disordered ground states for the two-circle potential defined by Eq. (26), 
varying the dimensionless parameter, a, while keeping fixed x = 0.4461. Large void 
regions emerge for the range a > 1. The void properties noticeably change as this 
parameter is increased. 



is a minimal energy configuration of the potential in two dimensions. Within 
each chain, however, the number of particles and relative order (in one dimen- 
sion) vary, with the particles being subject to an effective one-dimensional 
potential along the fcy-direction with an effective x value. This effective po- 
tential arises as a result of the structure factor being minimized on a set of 
unconstrained wave vectors (i.e., the parts of the black regions in the struc- 
ture factor plot of Fig. 3c which do not include the two circles). Taking a slice 
along a fixed ky (for all ky) gives the structure factor for a one-dimensional 
ground state subject to a potential which is a "donut" in one-dimensional 
/c-space, which is the aforementioned effective potential. The two topmost 
chains in Fig. 3c, being ordered, have high effective x values and the bottom 
two chains, being disordered, have low effective x values in this reduced di- 
mensional space. The variation comes as a result of the fact that the union 
of any two ground-state configurations is itself a ground-state configuration 
and that these one-dimensional structures are similarly entropically favored. 
Since the system is primarily constrained along the x-direction, the chains 
are free to be placed essentially randomly along the y-axis. 
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a = 0.4849 a = 0.9698 a = 1.4548 




XXX 



a = 1.9397 a = 2.9095 a = 3.8794 




Fig. 6 Two-dimensional pair correlation function for the two-circle potential Eq. 
(26) plotted in the x-y plane as a function of dimensionless parameter a,x = 0.4461. 
The Q values here correspond to the ones indicated in Fig. 5. Black regions indicate 
anti-correlation, g2 < 1. Lighter regions indicate spatial correlation, g2 > 1. Each 
plot has been averaged over twenty configurations. (Color figure online) 



As the system collapses to these chain-like structures oriented along the 
x-direction, a pair of wedge-shaped regions along the collapsed direction in 
Fourier space (the j/-direction) is forced to have positive S'(k); see Fig. 3c. 
The angle of the wedge near the origin determines the nature of the fila- 
mentary structures. The larger the wedge angle, the more tortuous are the 
filamentary structures. As the wedge angle decreases, the filamentary struc- 
tures are forced to become flat. The presence of these wedges explains why 
these configurations do not arise for isotropic potentials such as the stealth 
potential flO] . 



6.3 Three classes of structures 



6.3.1 Ordered chains 



For one such class of chains, those which display crystallinity in the direction 
of the chain, we are able to analyze the problem exactly. We derive an ex- 
pression for the structure factor of an ordered chain containing m particles, 
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Fig. 7 Two-dimensional static structure factor plotted logarithmically for the two- 
circle potential Eq. (26) as a function of dimensionless parameter a with fixed 
X = 0.4461. The a values here correspond to the ones indicated in Fig. 5. Black 
regions indicate that ^(k) has been constrained to zero (to within numerical error). 
The lightest regions indicate that 5(k) has reached a maximum, S{k) ~ A*', where A'^ 
is the number of particles. Each plot has been averaged over twenty configurations. 
(Color figure online) 



1e-10 



parallel to the a;-axis: 

\m\' 



m 
Ak 



e ) ( Z^e 1 

n=l ^ ^ n=l ^ 



cos{kxXom) — 1 
cos{kxXo) — 1 
cos(27rja;) - 1 

cos(27rj2;/m) - 1" 



(27) 



Therefore, there will be Bragg peaks at wave vectors nm/{2TTL), where n is 
an integer, m is the number of particles in the chain and L is the length of 
the fundamental cell along the direction of the chain. 

The form of the potential sets a bound on the number of particles which 
can form such a structure since no Bragg peaks can be found within the 
exclusion region. We call the fca;-distance from the origin to the first Bragg 
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peak I, the A:a;-distance from the origin to the perimeter of the exclusion 
region / and the diameter of the circle d. There are then two regimes in 
which an ordered chain can exist: 



f>l> 



l> f + d 
f + d 



(28) 
(29) 



The parameters d and / can be related to the parameters of the model, a 
and Xj and / can be related to the number of particles in the chain, giving 
us: 



a 



m > 



TT 1 
— > TO > - 

X 2 



— (a 

X 

— (a 

X 



2) 



2). 



(30) 



(31) 



A plot of an ordered chain and its structure factor is presented in Fig. 8, 
illustrating that the configuration minimizes the structure factor on wave 
vectors which fall within the two-circle exclusion region and is therefore a 
valid ground state. 




1eO 



1e-5 



'le-10 



Fig. 8 Sample ordered subconfiguration with A'^ = 31 particles (left panel) and 
its two-dimensional structure factor S'(k) plotted logarithmically. Black regions 
indicate that ^(k) has been minimized. The structure minimizes S(k) for the wave 
vectors within the non-zero region of the two-circle potential, indicating that the 
configuration is a ground state. The thickness of the vertical lines present in the 
structure factor is notable, indicating that the chain is not perfectly aligned in the 
j/-direction. (Color figure online) 



6.3.2 Ordered sinusoids 

We numerically calculate the structure factor for a representative sinusoid. 
For a sinusoid to be a valid ground-state structure, its structure factor must 
be minimized on the set of wave vectors specified by the potential. A plot of 
the numerically calculated structure factor is presented in Fig. 9. 
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Fig. 9 As in Fig. 8 except for a sample sinusoidal subconfiguration witii N = 87 
particles. Note that the period of the sinusoidal chain is twice the box length. (Color 
figure onhne) 



6.3.3 Disordered chains 

Analytical expressions for the structure factor of the disordered chains were 
difficult to derive, so numerical methods were employed again. A plot of the 
numerically calculated structure factor is presented in Fig. 10. 




Fig. 10 As in Figs. 8 and 9 except for a sample disordered configuration with 
N = 220 particles. (Color figure online) 



7 Concluding Remarks 

The ground-state properties for classical many-particle systems subject to 
a directional pair interaction were investigated. We found approximate re- 
lations for the pair statistics of the disordered ground states, which display 
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large voids with a tunable structure. The most ordered configurations display 
a collapse of dimensionality with the particles forming a few parallel, approx- 
imately linear, sometimes periodic chains in the x-direction with aperiodic 
placement of chains in the y-direction. This behavior is new and completely 
unexpected given previous studies of the collective-coordinate potentials. It 
should be noted, however, that so-called "patchy particles," with directional 
real-space interactions can also exhibit a collapse of dimensionality P^20j . 
However, such particles possess internal rotational degrees of freedom that 
are absent in the many-particle systems that we consider in this paper. 

We have explored the relations between isotropic and anisotropic pair po- 
tentials and their associated classical ground states, both here and in previous 
papers, highlighting the diversity produced by the collective-density variable 
approach [^[MB HTUlfTT] . However, many future avenues for research remain 
within these topics. As previously noted, each interaction potential generates 
a highly degenerate classical ground state occupying a domain whose dimen- 
sionality depends on the number of collective coordinate constraints that 
have been imposed. It will be important to devise mathematically precise 
descriptions of the connectivities, local curvatures, and configuration-space 
dispersions of these domains. Furthermore, extension of the present study 
of the two-circle potential to a three-dimensional analog might shed light 
on the collapse of particle arrangement dimensionality in the plane described 
above. Indeed, preliminary results for the corresponding "two-sphere" poten- 
tial model in three dimensions show once again that the particles exhibit a 
strong tendency to organize themselves into directed filamentary structures. 



A AfRne Theorem and the Ellipse Potential 

We derive a simple relation transform describing how the Fourier transform of 
a function in d-dimensional Euclidean space R"* behaves under invertible afHne 
transformations. We also present selected results for the ellipse potential described 
in Sec. III. The ellipse potential is described by the dimensionless aspect ratio, b/a, 
where b is the length of the minor axis and a is the length of the major axis of the 
potential in k-space. We go on to show that the ellipse potential is related to the 
previously studied circle potential by a simple afHne transformation, which allows 
for a simple analysis of the ground state patterns. 



A.l Affine Theorem for the Fourier Transform 

We derive a simple relation transform describing how the Fourier transform behaves 
under invertible affine transformations, generalizing a result in two-dimensions to 
arbitrary dimension [28J. Given a square-integrable function /(r) in d-dimensional 
Euclidean space, we can define its Fourier transform as: 




(32) 



The inverse Fourier transform is then naturally defined to be: 



/(r)=^r'[F(k)] 



(33) 
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Let /(r) be a function, / : E'' — !• R with a known Fourier transform T[f{r)] — /(k). 
Let us apply an invertible linear transformation, A, to r G R'', and denote the 
inverse by A~^. In order to be completely general, we also apply a translation 
denoted by the addition of a vector b. We now have: 



T[f{Av + h)]= I dr e-*"7(^r + b), 



where we have switched to linear algebra notation for ease and clarity of presen- 
tation. Let r' = Ar. In transforming the equation into the new coordinate system. 



we gain a factor 



from the Jacobian. 



Let k' = (k^A' 



dct(A) 

T[f{Ar + b) 

^f = {A-Y'^- 



yjdr' e---V;(r'+b). 



det(A) 



det(A) 
1 

det(A)' 
1 

det(A)' 
1 

det(A)' 



dr' e*'^"''/(r' +b) 



dr' e 



fir') 



'/(k') 



/((^"')'k), 



(34) 



where the second step follows from the Fourier shift property [25]. To the best of 
our knowledge, this result is the first of its kind. It is consistent with the previously 
mentioned two-dimensional result, as well as the standard results corresponding to 
multiplication by a constant factor and rotations '25 . 

Since our simulations are carried out in two dimensions, we will explicate the 
details of the previous result in this setting. For some invertible affine transforma- 
tion in two dimensions: 

r + To, 

where ^4 is a linear transformation given by the invertible matrix: 



A = 



a b 
c d 



and To is a translation vector: 



the Fourier transform of f[Ar) is: 



1 Zk-Tr 

a' 



(35) 



where /(k) is the Fourier transform of /(r). A' is a linear transformation defined 
by the matrix: 

f A ^.^1 (36) 



and A — det(yl) — ad — be. 
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In the case of the eUipse potential, the region on which C(k) is constrained can 
be expressed as a circular exclusion region to which we have applied invertible affine 
transformations. Since the ground-state properties of the circular exclusion region 
are well known (in light of previous studies), we can use this information, combined 
with this affine relation to draw conclusions about the ground-state properties of 
the ellipse potential. We emphasize that the affine transformation is not applied to 
the particles themselves and therefore does not allow for the particles to rotate. 




Fig. 11 Crystalline ground state configuration for the ellipse potential. If the con- 
figuration is subject to the squeeze transformation described above, the triangular 
lattice is not recovered. 



A.2 Results 

A. 2.1 Disordered ground states 

Ground state configurations were generated over a range of aspect ratios for the 
ellipse potential at a fixed x value. Comparing the pair correlation function of 
the ground states of a potential of given aspect ratio to that of the isotropic circle 
potential, it is apparent that the approximate affine relation derived in the previous 
section holds. In addition, if x is increased to values near but still below 0.5, the 
particles develop what resemble hard cores and we see the emergence of perfect 
nematic order (which is by construction). 



A. 2. 2 Crystalline ground states 

One might hypothesize that the crystalline ground states would obey the same affine 
relationship that the potential exactly obeys and the pair statistics approximately 
obey, i.e., that the absolute ground state would be an appropriately transformed 
triangular lattice. However, this is not the case. As can be visually ascertained from 
Fig. 11, applying the inverse squeeze transformation does not recover the triangular 
lattice. In order to confirm this, a perturbed transformed lattice configuration was 
used as an initial particle configuration for simulations using the ellipse potential. 
Not only did the trials not find the desired lattice to be the ground state, but no 
true ground state was reached when these initial conditions were used. However, it 
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should be noted that 418 particles cannot be used to created a nearly undeformed, 
affine transformed triangular lattice in a square box (as is possible with the un- 
transformed triangular lattice), which may have influenced the convergence of the 
procedure. 
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